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SECTION  I 

INTRODUCTION  AND  SUMMARY 


The  motivation  for  the  present  investigation  was  to  see 
if  reorienting  the  crystal  direction  of  silicon  photomulti- 
plying  detectors  could  lead  to  a  reduction  in  their  noise  at 
high  gain.  The  mean  square  current  fluctuations  in  such  de¬ 
tectors  is  proportional  to  the  gain  squared  as  long  as  the 
gain  is  less  than  the  ionization  rate  ratio.  If  the  gain 
exceeds  the  ionization  rate  ratio  then  the  noise  increases 
to  something  like  the  third  power  of  the  gain.  Of  all  the 
semiconductors  investigated  so  far,  silicon  has  the  highest 
ionization  rate  ratio,  ranging  from  five  to  one  hundred  de¬ 
pending  on  the  value  of  the  electric  field  in  the  avalanche 
region. 

Calculations  of  the  threshold  energies  in  various  crys¬ 
tal  directions  had  been  published  by  Anderson  and  Crowell1 
showing  nearly  equal  thresholds  of  about  3  eV  for  both  elec¬ 
trons  and  holes  for  the  <Clll>  direction.  This  is  the  direction 
for  which  the  high  ionization  rate  ratio  is  observed.  Fur¬ 
ther,  their  calculations  showed  that  in  the  CLOCi)  direction 
the  threshold  voltage  for  holes  could  be  twice  that  ( to  1  eV ) 
of  tne  electrons.  In  tne  Ol0>  direction  the  ionization 
thresholds  were  about  the  same  at  approximately  2  eV,  so 
presumably,  the  rates  in  this  direction  would  be  similar  to 
the  ratio  found  in  the  (ill)  direction. 

The  results  of  the  measurements  reported  here,  contrary 
to  exoectations,  show  the  highest  rates  and  the  smallest 
ratio  of  rates  is  obtained  in  the  <(l00)>  crystal  direction, 
a  slightly  larger  ratio  is  obtained  in  the  <(llO>  direction, 
and  the  largest  ratio  is  obtained  in  tne  <111}  direction. 

Since  a  careful  recnecking  of  the  measurements  yielded  the 
seme  result,  some  thought  has  been  given  to  the  threshold 


calculations.  In  Anderson  and  Crowell's  work  only  colmear 
momentum  constructions  of  the  initial  and  final  particles 
through  the  zone  center  (i.e.  k=0)  were  considered.  During 
the  course  of  this  work  colinear  momentum  constructions  going 
througn  the  conduction  band  minima  and  not  through  the  zone 
center  were  considered.  Neither  of  tnese  sets  of  construc¬ 
tions  corresponds  to  the  observed  ionization  rates.  Non- 
colmear  constructions  have  been  considered  here  in  which  the 
final  particle  velocities  are  colinear,  but  the  wave  vector 
from  the  Lrillouin  zone  center  to  the  conduction  band  mini¬ 
mum  and  the  initial  wave  vector  need  not  be  colinear  with  the 
final  velocities.  This  non-colinear  problem  is  much  more 
complex  and  is  not  solved  in  tnis  report,  but  an  approximation 
of  it  is  presented  that  qualitatively  agrees  with  the  observed 
orientation  dependence. 


SECTION  II 

DESCRIPTION  OF  EXPERIMENT 


1.  Sample  Preparation 

A  p-type  boule  of  silicon  was  oriented  and  wafers  were 
cut  in  the  three  principv  .  crystallographic  directions  <100>, 
<110),  and  <111> .  The  measured  resistivity  was  .49.ft.-cm  for 
the  <100>  sample,  .68,ft-cm  for  the  <110>  sample,  and  .50&-cm 
for  the  <111>  sample.  The  samples  were  then  diffused  at 
about  HOC  °C  in  an  induction  heated  vacuum  furnace  for  lb 
to  60  minutes.  The  tantalum  bucket  enclosing  the  wafers  is 
prebaked  at  about  1700  °C  to  get  rid  of  electrically  active 
impurities.  A  diffusion  grade  phosphorous  nitride  was  used 
as  a  donor  source.  The  contamination  background  of  this 
diffusion  system  is  less  than  10  /cm  added  impurities. 

The  time  and  the  temperature  were  chosen  to  produce  linearly 
graded  junction  impurity  profiles.  After  diffusion  some  of 
these  wafers  were  thinned  to  a  thickness  of  about  one  mil  so 
that  hole  and  electron  injection  in  the  junction  could  be 
accomplished  from  opposite  sides  of  the  wafer. 

2.  Diode  Profiling 

The  electric  field  profile  of  the  junctions  is  obtained 
from  capacitance  measurements  as  a  function  of  reverse  bias. 

The  data  for  a  <110>  diode  is  shown  in  Fig.  1  where  it  is  seen 
that  the  inverse  cube  of  the  capacitance  is  a  linear  function 
of  the  bias  out  to  breakdown.  This  linearly  graded  profile 
allows  a  more  accurate  determination  of  the  maximum  field  in 
the  junction  than  an  asymmetric  step  junction.  The  reason 
is  tnat  in  tne  step  junction  free  carrier  spillage  from  the 
heavily  doped  side  significantly  affects  the  maximum  junction 
field.  Thus  for  a  linearly  graded  junction  the  simple  Schottky 
approximation  can  be  used  to  accurately  determine  the  electric 
field.  The  following  constants  were  determined  for  the  junc¬ 
tion  of  Fig.  1: 
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Figure  1.  Diode  capacitance  data  in  the  <liO>  direction. 


Area  =  4.2  x  10  ^  cm2  , 

dN  ,  ,,.21  -4  , 

—  =  1.73  x  10  cm  ’ 

V,  .  =  0.7  Volts 
bi 

=  4.2  x  104  (V  +  Vb.)2/3  .  (1) 

The  capacitance  data  for  the  <100>  oriented  diode  is 
shown  in  Fig.  2.  This  data  shows  some  curvature  which  is  a 
consequence  of  the  fact  that  the  diffusion  time  was  not  long 
enough.  The  phosphorous  distribution  has  roughly  an  expon¬ 
ential  distribution  into  a  uniform  background  of  about 
16  3 

4.0  x  10  /cur  acceptor  concentration.  The  electric  field 
near  the  maximum  still  has  a  parabolic  shape,  but  the  electric 
field  dependence  on  the  voltage  is  altered  significantly  near 
breakdown.  The  solid  line  in  Fig.  2  shows  the  fit  obtained 
by  assuming  an  exponential  diffused  impurity  gradient  and 
a  uniform  background  of  acceptors, 


N(x)  =  Nq ( exp(-x/L)  -  1)  .  (2) 

Using  this  impurity  distribution  one  can  derive  the  follow¬ 
ing  expressions  for  the  electric  field  and  the  total  junction 
potential : 
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where  u>  =  w/L,  and  w  is  the  junction  depletion  width.  The 
parameters  of  the  fit  to  the  data  shown  in  Fig.  2  are: 

_3  2 

Area  *  1.7  x  10  cm"  , 

„  .  ,  iri  16  4 

N  =  4.1  x  10  cm  , 

o 

L  =  1.28  x  10  cm  , 

-  0.6  Volts 


Note  that  the  background  doping  agrees  well  with  that  deter- 
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Diode  cap.-iei t.nnce  data  in  the  <100)  direction 


mined  from  resistivity  measurements  before  diffusion. 


Capacity  data  for  the  <111>  direction  is  not  presented 
because  multiplication  data  were  directly  compared  with  cal¬ 
culated  multiplication  using  previously  published  ionization 

p 

rates.  In  addition,  all  of  the  parameters  given  above  for  the 
<100>  and  the  <110)  have  an  error  tolerance  of  about 

3.  Photomultiplication  Measurements 

The  electron  and  hole  multiplication  data  for  the  <110) 
and  the  <100>  directions,  which  are  shown  in  Fig.  3,  were  ob¬ 
tained  using  two  different  techniques.  To  obtain  pure  elec¬ 
tron  and  hole  excitation  of  the  avalanche  a  <110)  oriented 
wafer  was  thinned  to  a  thickness  of  slightly  less  than  25  /im, 
and  non-penetrating  radiation  was  used  to  illuminate  the  diode 
from  opposite  sides  of  the  junction.  The  multiplication  data 
obtained  this  way  was  compared  with  that  obtained  by  using  non¬ 
penetrating  radiation  (435.8  nm  interference  filter)  from  the 
n-side  for  hole  excitation  and  using  a  GaAs  filter  to  obtain 
penetrating  radiation,  again  from  the  n-side,  to  obtain  the 
electron  excitation.  This  second  technique  relies  on  the  fact 
that  only  a  small  number  of  holes  will  be  generated  in  the  dif¬ 
fused  layer  of  about  1  micron  thickness  while  electrons  are  coll¬ 
ected  over  a  distance  of  some  30  microns  at  least  in  the  p- 
type  substrate  when  penetrating  radiation  is  used.  Both 
techniques  gave  multiplication  data  for  electrons  and  holes 
that  agreed  to  within  about  5c/°,  so  the  simpler  second  technique 
was  used  to  obtain  the  data  on  the  <100>  sample. 

For  comparison,  Fig.  4  shows  calculated  electron  and  hole 
multiplication  in  the  <111)  direction  for  a  breakdown  voltage 
close  to  that  of  the  <110)  data.  Perhaps  the  best  way  to 
illustrate  the  differences  between  the  three  crystal  direc¬ 
tions  is  to  plot  the  multiplication  versus  the  normalized 
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Figure  i.  Multiplication  data  for  <110)  and  <100)  directions. 


8 


voltage  as  is  shown  in  Fig.  5.  Here  tne  highest  multiplica¬ 
tion  and  the  smallest  ratio  of  electron  and  nole  multiplication 
is  obtained  in  the  <100>  direction.  In  the  <110>  direction 
the  multiplication  is  somewhat  less  and  the  ratio  of  electron 
to  hole  multiplication  is  slightly  larger.  In  the  <111>  dir¬ 
ection  the  electron  multiplication  is  seen  to  o<j  slightly 
smaller  than  the  hole  multiplication  in  the  <100>  direction 
wmle  the  hole  multiplication  is  very  much  lower.  This  data 
implies  that  the  lowest  threshold  ionization  energy  should  be 
founa  in  tne  <100)  direction,  an  intermediate  threshold  for 
the  < 110)  direction  while  the  highest  tnreshold  should  be 
found  for  the  <111>  direction.  The  implication  of  this  data 
is  thus  at  variance  with  the  thresholds  calculated  by  Ander¬ 
son  and  Crowell  except  for  the  electron  threshold  m  the  <100> 
direction.  Ine  last  section  will  consider  non-colmear  mom¬ 
entum  constructions  that  can  account  for  the  observed  multi- 


Figure  5.  Multiplication  data  normalized  to  breakdown  voltage. 
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SECTION  III 
IONIZATION  RATES 


For  the  linearly  graded  junctions  we  are  considering  the 

integral  equations  relating  tne  ionization  rates  to  the  muiti- 

2 

plication  data  are  given  by 
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These  integrals  are  special  cases 
and  the  solutions  to  Eqs.  (5)  and 
following  form  which  is  suitable 
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3 

of  Abel's  integral  equation 
(6)  can  be  written  in  the 
for  numerical  evaluation: 
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/3<V 


ds( l-s2)T 


(8) 


where  6=  6  (1  -  s2) ,  and  the  voltage  derivative  of  the 
m 

function  of  the  multiplication  is  evaluated  at  the  voltage 

Vg  corresponding  to  a  maximum  electric  field  m  the  junction 

of  i  -  6(1  -  s2). 
s  m 


The  multiplication  data  shown  in  Fig.  3  has  been  used 
to  evaluate  Eqs.  (7)  and  (8).  It  has  been  noted  that  the 
profile  of  the  <100)  sample  is  not  quite  a  linearly  graded 
junction,  but  it  has  been  estimated  that  the  slight  asymmetry 
of  the  electric  field  distribution  in  the  region  where  the 
ionization  rates  are  appreciable  should  not  change  the  result 


1? 


obtained  by  assuming  a  linear  gradient  by  more  than  The 

calculated  ionization  rates  are  plotted  in  Fig.  6.  For 

comparison  the  previously  published  rates  for  the  (111)  dir- 
2 

ection  are  shown  in  Fig.  7.  At  the  present  time  it  is  not 
clear  why  the  slope  of  the  ionization  rates  gets  progress¬ 
ively  more  negative  as  one  goes  from  the  (100)  to  (110)  to 
the  (111)  direction. 
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Figure  6.  Ionization  rates  for  <100>  and  <I10>  airections 
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SECTION  IV 


IONIZATION  THRESHOLD  ENERGIES 
1.  Old  Threshold  Estimates 

The  most  complete  estimate  of  the  ionization  threshold 
energies  was  carried  out  by  Anderson  and  CrowelIx  wr.o  pro¬ 
ceeded  from  the  pseudopotential  calculations  of  Chelikowsky 
and  Cohen,  which  are  shown  in  Fig.  8  along  with  the  zone 
scheme  for  silicon  which  is  shown  in  Fig.  9  for  reference. 

In  working  out  their  graphical  solutions  they  assumed  that 
the  wave  vector  of  the  initiating  particle,  the  velocities  of 
the  three  particles  after  the  collision,  and  the  wave  vector 
from  the  zone  center  to  the  conduction  band  minimum  were  all 
colinear.  Figure  ID  shows  the  <100>  trajectory  (X  -  T  -  X), 
and  the  <110> trajectory  (X  -  K  -  F  -  K  -  X) .  Figure  11 
shows  the  trajectory  in  the  <111>  direction  (L  -  T  -  L) . 

Their  lowest  threshold  energy  is  for  electrons  in  the  <100> 
direction,  1.1  eV .  The  threshold  energy  for  holes  in  the  same 
direction  was  1.8  eV .  In  the  <110>  direction  they  found  a 
threshold  energy  of  3*1  eV  for  electrons  and  a  hole  energy 
of  1.8  eV.  Finally  in  the  <111>  direction  they  obtained 
threshold  energies  of  3*1  eV  for  electrons  and  2.9  eV  for 
holes.  As  noted  before  these  threshold  energies  do  not 
correspond  to  the  observed  multiplication. 

In  the  course  of  this  work  it  was  realized  that  other 
trajectories  should  be  considered.  Additional  trajectories 
were  considered  that  were  parallel  to  the  principal  crystal 
directions,  but  passed  through  the  conduction  band  minimum 
along  the  <100)  axes.  The  band  energies  along  these  trajec¬ 
tories  were  plotted  out  out  the  estimated  threshold  energies 
did  not  even  qualitatively  agree  with  the  observed  data. 


ie.-’ 


Figure  10.  Particle  trajectories  through  Brillouin  zone  in 
<110>  and  <100>  directions  passing  through  the  zone 
center  and  through  the  conduction  band  minimum. 
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Figure  II.  Particle  trajectories  through  the  Bn^Iouin 
20.it;  in  the  <L11>  direction  passing  through  z:\e  zone 
oeni/cr  and  through  the  conduction  bond  minimum. 


These  trajectories,  however,  also  made  the  assumption  tnat  all 
the  vectors  were  colinear.  It  was  eventually  realized  tnat 
the  wave  vector  of  the  initial  particle  need  not  be  colinear 
with  the  final  wave  vectors  nor  with  the  wave  vector  from  tne 
center  of  the  zone  to  the  conduction  band  minimum.  Since  the 
geometry  can  be  quite  complex,  consider  the  situation  illus¬ 
trated  in  Fig.  12.  Here  the  k^-axis  is  coincident  with  the 
<100>  direction  and  the  k  -axis  intersects  the  <111>  direction. 

y 

We  also  consider  an  electric  field  parallel  to  the  <111>  dir¬ 
ection,  making  the  angle  psi  equal  to  54.7  degrees.  The  init¬ 
ial  velocity  will  not  be  parallel  to  the  field,  but  will  make 
an  angle  6  with  the  ky-axis.  Since  the  secondary  hole  is  in 
a  nearly  spherical  band  its  velocity  and  wave  vector  are  essen¬ 
tially  parallel.  Then  we  note  that  the  extremum  of  minimum 
threshold  energy  unaer  the  condition  of  conservation  of  crys¬ 
tal  momentum  requires  that  all  three  final  velocities  ne  co- 
lmear.  Thus  if  the  hole  is  in  a  spherical  band  while  the  two 
electrons  are  in  a  band  with  ellipsoidal  energy  surfaces, 
the  wave  vectors  will  not  in  general  be  parallel. 

2.  New  Processes  for  Determining  Threshold  Energies 

We  have  posed  a  difficult  problem  which  has  not  been 
solved,  but  there  is  an  approximate  solution  which  qualita¬ 
tively  agrees  with  the  observed  multiplication  data.  The 
approximation  we  wish  to  consider  is  the  following:  at  the 
minimum  threshold  energy  the  particle  energies  after  the  colli¬ 
sion  should  be  quite  small,  thus  only  a  small  error  should  be 
incurred  if  one  assumes  the  wave  vectors  and  the  velocities  of 
the  final  particles  to  be  parallel.  The  wave  vector  from  the 
zone  center  to  the  conduction  band  minimum,  however,  is  large 
and  we  still  allow  this  vector  to  be  non-colinear .  We  also 
will  assume  that  the  band  energies  are  parabolic  in  the  plane 
containing  the  electric  field  and  passing  through  the 


conduction  band  minimum  or  the  valence  band  maximum. 


The  assumption  of  parabolic  bands  allows  us  to  write  the 


energies  of 

the  initial  and  final 

three  particles, 

E2 

,=  *2(k2,-  ko)2/2me 

1 

(9) 

E1 

=  *2(V  ?o)2/2me 

(10) 

E2 

=  /h2(k2-  ko)2/2me 

» 

(11) 

E1 

,=  /h2(k1,)2/2mh 

* 

(12) 

where  the  subscript  2'  denotes  the  initial  electron  (in  this 
case),  1  and  2  denote  the  secondary  electron  and  the  initial 
electron  after  collision,  and  1*  denotes  the  secondary  hole. 
Conservation  of  energy  and  crystal  momentum  give  the  follow¬ 
ing  relations, 

E2,=  Ex.  E2*  Elt.  Eg  ,  (15) 
*2'=  *1+  *2+  *1’  *  (14) 


where  E  denotes  the  energy  gap  and  phonon  cooperation  is 
S 

neglected.  The  minimum  threshold  energy  is  obtained  from 
the  variational  equation 


(15) 


where  the  vector  0  is  a  Lagrange  multiplier.  The  result  of 
Ea.  (15)  is  that  all  the  final  wave  vectors  are  colinear, 


k 


1 


k 


o 


(me/mh> 


(16) 


and  in  the  more  general  case  it  would  be  the  velocities  which 
would  be  colinear.  The  non-colinear  part  of  the  problem,  ky, 
not  parallel,  in  general,  to  the  wave  vectors  k9,  ,  k^,  k„,, 
or  k^, .  By  a  somewhat  lengthy  manipulation  of  Eqs.  (13), 


is 
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( 14 )  and  (15)  we  can  obtain  an  expression  for  E.?,  in  terms  of 
the  band  gap  energy,  the  displacement  of  the  conduction  band 
minimum  from  the  zone  center,  and  the  ratio  of  the  hole  to 
electron  effective  mass, 


-->  i 


;  2+  H- 

~  1+fi 


1+  -*-2cos  6 


(1+/*  )* 


2cos6 


jx^x(l+fi+cos^©)  +  (2i-fi  )  (1+  fj.  )j  |  ,  (17) 


(1+  H- ) 

where  9  is  the  angle  between  the  wave  vectors  k 


L  and  k  , 

^  p  J  2 

ll  is  tne  mass  ratio  m./m  ,  and  x  =  E  /E  ,  where  E  =  A  k  /2 
“  h  e  o  g  oo 


m 


g  o  •  ' o  '  e 

It  is  immediately  apparent  that  Eq.  (17)  reduces  to  a  familiar 


result  for  a  direct  gap  semiconductor  (i.e.  kQ=  0) 


'2' 


2^ 
1+  ft 


(18) 


Another  case  of  immediate  interest  is  the  < 100>  direction  in 
silicon  where  9  =  0  and  all  the  vectors  are  colmear 


2±JL 


2±JL 


E 


(1+fO 


2 


x  - 


(1+  ft  ) 


2  L 


x(x+l+  ft  )  (2+/i  )] 


(19) 


Equation  (19;  contains  the  surprising  result  that  for  x  =  1 
trie  threshold  energy  is  exactly  the  band  gap  energy  whereas 
for  x  =»  0  we  get  the  classical  result  of  Eq.  (18).  One 
would  expect  that  displacing  the  conduction  bana  minimum 
away  from  the  zone  center  would  increase  rather  than  decrease 
the  tnreshold  energy.  One  further  point  is  to  be  made  concern¬ 
ing  Eq.(17)  and  that  is  that  one  can  obtain  the  case  of  a 
hole  being  the  energetic  initial  particle  by  simply  invert¬ 
ing  the  mass  ratio  ^  .  This  fact  is  made  more  transparent  by 
translating  the  origin  in  k-space  to  the  conduction  bana 
minimum  when  an  energetic  hole  is  the  initial  particle. 
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We  now  apply  Eq.  (17)  to  the  three  crystalline  directions 
in  silicon  by  calculating  the  effective  mass  m  a  given 
orientation  with  respect  to  the  energy  ellipsoids  from  the 


equation 


f?  k  2 

_ x 


b\2 

2m  * 
e 


and  the  value  of  mg*  is  also  used  to  evaluate  the  expression 

for  x  o  q 

£  h  ^rf  n 

x  =  —  =  - —  •  —  (21) 

6  e  g 

and  where  kQ  is  taken  to  be  0.9  times  the  distance  to  the 

band  edge  in  the  <100>  direction.  Using  the  effective  masses 

m  =  (.98)m  ,  m.*  (.19)m  .  and  m  =  (,92)m  we  have  calculated 
1  o  t  o  h  o 

the  threshold  energies  of  electrons  and  holes  for  the  three 
directions  and  listed  them  in  Table  1. 

It  is  clear  from  the  calculated  thresholds  listed  in 
Table  I  for  the  three  crystal  orientations  that  they  are 
consistent  with  the  experimental  measurements.  In  fact  the 
agreement  is  sufficiently  encouraging  to  warrant  carrying  out 
the  complete  calculation  where  the  velocities  and  wave  vectors 
are  not  assumed  to  be  colinear. 


TABLE  I 


CALCULATED  THRESHOLD  ENERGIES  IN  SILICON 


Orien. 

0 

degrees 

X 

H- 

E2,/-g 
elec . 

E2'/ES 

holes 

<  100) 

0. 

.93 

.93 

1.00 

o 

o 

90. 

9.03 

2.74 

2.6; 

- 

9-42  j 

<110) 

49. 

3.01 

1.63 

i.  41 

1.96 

90. 

9.03 

2.74 

2.62 

9.42 

<111> 

94.74 

3.68 

2.00 

1.62 

_ 

2.02 
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